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Abstract 

We present the details of a global analysis of electroweak data, including fermion 
masses and mixing angles, in SO (10) SUSY GUTs. Just as precision electroweak data is 
used to test the Standard Model, the well determined Standard Model parameters are 
the precision electroweak data for testing theories beyond the Standard Model. In this 
paper we use the latest experimentally measured values for these parameters. We study 
several models discussed in the literature. One of these models provides an excellent 
fit to the low energy data with ~ 1 for 3 degrees of freedom. We present graphs of 
constant contours as functions of position in soft SUSY breaking parameter space, 
as well as our predictions for a few selected points in parameter space. We also study 
the sensitivity of our results to changes in various parameters. Finally, we discuss the 
consequences of our work in the context of a general MSSM analysis at the Z scale. 
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1 Introduction 



The Standard Model is by far one of the greatest achievements of the twentieth century. Its 
stature is further elevated by the fact that after a decade or more of testing there is still no 
evidence for physics beyond the Standard Model(SM). Yet we nevertheless remain firm in our 
belief that the SM cannot be the whole story. For one thing there are 18 phenomenological 
parameters which await derivation from a more fundamental theory. In particular 13 of 
these parameters are related to the fermion mass sector. Secondly, there is growing evidence 
that neutrinos may have mass; yet another indication of physics beyond the SM. Finally, 
a quantum theory of gravity has yet to be incorporated into the SM. Clearly string theory 
is the leading candidate for a fundamental theory of Nature; it can incorporate gravity, as 
well as the known quarks, leptons and gauge interactions. The problem lies in finding the 
matching conditions between this fundamental string theory and the effective field theory 
which describes the physics below the Planck (or string) scale.Q 

The fermion mass sector is now well measured. It is clear that fermions come in well 
defined families with an amazingly simple regularity in the hierarchy of masses and mixing 
angles. It is hoped that by understanding this sector we can find the correct matching 
conditions at the string scale. Now seems an appropriate time to begin such an endeavor. It 
was only ten years ago that the top quark mass and the CKM angle Vub were unknown; while 
Vcb and the weak mixing angle sirt^Ow were measured, but with large error bars. Today, as a 
result of both experimental and theoretical progress, all but one of the 18 parameters of the 
Standard Model are known to much better accuracy [the Higgs mass has yet to be measured]. 
These parameters are the precision electroweak data for testing theories of fermion masses. 

Of course once setting out to test theories of fermion masses, we must necessarily choose 
some particular theories to test. It is in this choice of a theory that we now invoke some 
theoretical prejudices. It seems clear to us that the simple observed pattern of fermion 
masses and mixing angles is not due to some random dynamics at an effective cut-off scale 
M but is instead evidence that a small set of fermion mass operators are dominant. Thus 
we are lead to postulate - 

1. that a few effective operators at an effective cut-off scale M (where M = Mpianck or M string) 
(or at a GUT scale Mq) dominate the quantitative behavior of fermion masses and 
mixing angles; and 

2. that a more fundamental theory, incorporating Planck or stringy dynamics, will gen- 
erate an effective field theory below M including these dominant operators. 

We are thus encouraged to find this effective field theory, thereby determining the match- 
ing conditions near the Planck scale. However, in order to make progress, it is clear that 

^Within string theory the problem is different; it is to understand string dynamics and why one particular 
vacuum is preferred over all others. 
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order of magnitude comparisons with the data are insufficient. Moreover, since the predic- 
tions of any simple theory are correlated, a significant test of any theory requires a global fit 
to all the low energy data. This paper makes the ffist attempt to bring the tests of theories 
of fermion masses into conformity with the accuracy of the low energy data. 

In this paper we present a global analysis of precision electroweak data, including 
fermion masses and mixing angles, within the context of several theories of fermion masses 
based on SO(IO) SUSY GUTs. Of course, the self-consistent global analysis we describe can 
be applied to any predictive theory. 

Why look at SO(IO) SUSY GUTs? 

• We use SUSY GUTs since they give the simplest explanation for the experimental 
observation that the three gauge couplings appear to meet at a scale of order 10^^ 
GeVi. 

• We use SO(IO) since it provides the simplest explanation for the observed family struc- 
ture of the light fermionsp]. 

An important feature of our analysis is the inclusion of low energy supersymmetric thresh- 
old corrections. We study their significance for minimizing x^- As has been shown earlier 
with tan/5 large there are potentially large corrections to those observables which are 
tan/3 suppressed at tree level. These observables include, for instance, masses of b, s and d 
quarks, and masses of charged leptons. The CKM elements Vcb and Vub, resulting from the 
mismatch between up- and down-type quark mass matrices, also get these corrections. In this 
case the large corrections to the d-quark mass matrix are passed on to the unitary matrices 
which diagonalize it 0]. All these corrections could potentially conspire and open a window 
in parameter space which has been unnoticed in the analysis neglecting them. In fact, such 
interplay is hard to study in the usual bottom-up scenario, but simple to implement in a 
global top-down analysis like this one. 

In section 2 we describe the low energy observables which enter the analysis, their 
measured values and the associated experimental and/or theoretical errors. 

In section 3 we discuss the self-consistent evaluation of the theoretical predictions, starting 
with the effective SO(IO) GUT and the resulting boundary conditions at Mq; then detailing 
the renormalization group running to Mz and the threshold corrections included at the weak 
scale. 

In section 4 we discuss models with four effective mass operators at the GUT scale. We 
show in detail results for the best model out of nine ADHRS|^ models. In this model [model 
4] the number of arbitrary Yukawa parameters riy = 5. The best fits give ~ 13 — 14 for 
5 degrees of freedom. All low energy observables are fit to within their 2a bounds. The fits 
for 4 observables, however, (Vcb, Vub/Vcb, the Kaplan-Manohar-Leutwyler ellipse parameter 
Q and the Jarlskog CP violating invariant J) lie close to the 2a border. We show that these 
results are true predictions of the theory, i.e. we see no other way to improve the results. 
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except with the addition of a new operator contributing to the 13 and 31 elements of the 
Yukawa matrices 0. 

In section 5 we analyze 3 models derived from complete SUSY GUTs discussed recently 
by Lucas and S.R.0]. One of these three models is identical to model 4 of ADHRS when 
states with mass greater than Mq are integrated out; this is model 4b. The other two include 
one new effective operator at the GUT scale which distinguishes the two models, 4a and 4c. 
We find that model 4a gives a best fit with ~ 3 for 3 dof in a small corner of parameter 
space, while model 4c gives <1 for 3 dof over a large subspace of the allowed parameter 
space. 

Our main results are found in figures [1 - 4] and tables (III - IX). These include graphs 
of constant contours in the two dimensional soft SUSY parameter space labelled by mg, a 
universal squark and slepton mass parameter, and M1/2, a universal gaugino mass parameter 
both with values given at Mq- We have several such plots for different values of the Higgs 
mass parameter /i, given at Mz- In the tables we give the computed values for the low 
energy observables for a few selected points in parameter space. We also include values for 
the Higgs and SUSY spectra, the CP violating angles a, j3 and 7, measurable in neutral B 
decays, the Wolfenstein parameters r], p, and the ratio mu/md for these points. Note, in 
figure 5, we present a scatter plot for values of (sin 2a, sin 2/?) predicted by models 4a and 
4c for points in figures 3 (a - c) and 4 (a - c), respectively, with values < 3. Finally we test 
the sensitivity of x^ to changes in various parameters. 

Section 6 contains further discussion of our results and conclusions. 



2 Low energy observables — experimental values 

Our x^ function includes 20 low energy observables. In addition we have incorporated the 
experimental bounds on sparticle masses into the code as a penalty in x^; added if one 
of these bounds is violated. This gaurantees that we remain in the experimentally allowed 
regions of parameter space. Let us now discuss the experimental observables and their errors. 
We have 

• 6 parameters associated with the Standard Model gauge and electroweak symmetry 
breaking sectors — as{Mz), aEM, G^, Mw, Mz, and pnew] 

• 13 parameters associated with fermion masses and mixing angles — Mt, mh{Mb), {M^— 
Mc), m^, md/rris, M^, M^, Mg, 14s, Kb, Kfe/Kb, and J ;0and 

• the branching ratio for b s'y, 

^Note that models with smaher values of tan/3 not only involve an extra free parameter but also require 
additional operators to correctly reproduce the experimental values for the first generation quark masses j^] . 

■^In the actual analysis J is replaced by the hadronic matrix element Bk as will be explained later in 
this section. 
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where pnew includes the added contribution to the electroweak p parameter from new physics 
beyond the SM, Q is the Kaplan-Manohar-Leutwyler elhpse parameter^] relating m, d, and 
s quark masses and J is the Jarlskog CP violating invariant (see more on these parameters 
below). Note also, that mass parameters denoted with a capital letter M are defined as pole 
masses, while m;,, mg, m^, mu are defined as the MS running masses. The three hght quark 
running masses are evaluated at the scale 1 GeV. 

The experimental values for the observables are given in table |1] with the associated 
experimental or theoretical error a. Note, a is taken to be either the experimental error or 
1/2%, whichever is larger. This is because our theoretical calculation also introduces an 
error; the numerical solution to the renormalization group equations is computed with some 
limited precision and the values of all observables are calculated within the first orders of the 
perturbative expansion in couplings. Thus 1/2% represents a conservative simple estimate 
of the combined theoretical error assigned to each low energy theoretical output As a 
consequence, the listed errors for Mz, Mw, G^, a em, Mr, M^, Me are dominated by 
theoretical uncertainties. 

For as{Mz) we use the central value quoted by Schmelling in Warsaw [p^]. However, for 
the error we choose a conservative value as suggested by Webber [|TT[ or Burrows |T^. This 



accounts for the fact that the systematic errors of individual as measurements are significant 
and thus caution is required when adding them together. 

Note that the traditional Standard Model parameter sin"^ 9 w{Mz)\jjs is not among 
the observables in table [^. Neither is J, ex or any other standard CP violating quantity 
listed. Instead, in the list of observables, sin^ 9w{Mz) is replaced by the Fermi constant 
Gfj_, and J is replaced by the bag constant I3k- is substituted for s\v? 6\y(Mz) since 
the former is extracted from the muon lifetime formula [|13| with very little sensitivity to 
the particle content at the scale Mz and the underlying theory at this scale. In fact, it 
could be defined completely in the context of the effective four-fermi contact interaction plus 
electromagnetism, and the remaining uncertainty of the order M^/M^ would then appear 
in the quantities derived from G^. On the contrary, the experimental value of siv? 6w{,Mz) 
in the MS scheme is derived from the precision measurements at the Z peak. Loop effects 
play an important role in its precise determination and hence one has to specify the theory 
at Mz and its particle content, and then calculate all relevant radiative corrections. Thus 
the underlying theory is an explicit input for the extraction of the experimental value of 
siv? 6w{,Mz) and, naturally, it is commonly taken to be the Standard Model (SM). However, 
we adopt an approach where the minimal supersymmetric extension of the Standard Model 
(MSSM) is matched directly to the low energy SU (3)c x U{l)em effective theory by integrating 
out the superpartners at the same time as the top quark and W/Z bosons (at a single scale 
Mz)', hence there is no room for the SM. This approach is particularly sensible if there are 
light SUSY particles (below Mz), such as charginos, neutralinos, pseudoscalar Higgs, etc.. 



"^For we take cr to be 1%. This is to account for the additional error resulting from our neglecting one 



loop SUSY box and vertex corrections. 
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for then there is no energy regime where the SM dominates |114[|. Although not among the 
observables contributing to the function, the SM value of sin^ 9wiMz) is displayed among 
the model predictions in tables ^ and |[ 

The CP violating parameter J is defined by the expression - 

J = Im{VudV:,VM) ^ iKdllW^^llK^r sine (1) 

where ^ is the CP violating phase. We test J by a comparison to the experimental value 
extracted from the well-known K*^ — mixing observable = (2.26 ± 0.02) x 10~^. The 
largest uncertainty in such a comparison, however, comes in the value of the QCD bag 
constant Bx- We thus exchange the Jarlskog parameter J for Bx in the list of low-energy 
data we are fitting. Our theoretical value of Bk is defined as that value needed to agree with 
ex for a set of fermion masses and mixing angles derived from the GUT-scale. We test this 
theoretical value against the "experimental" value of B^, which is given in table |l|. This 
value, together with its error estimate, is obtained from recent lattice calculations [0]. 

The experimental value for pnew is obtained from Langacker's combined fits to the preci- 
sion electroweak data, presented at recent workshops |1l6[ . 



For fermion masses and mixing angles we use those combinations of parameters which 
are known to have the least theoretical and/or experimental uncertainties. For example, 
while the bottom and top quark masses are known reasonably well, the charm quark mass 
is not known as accurately. On the other hand, heavy quark effective theory relates the 
mass difference Mh — Mc between the bottom and charm quark pole masses to about 5% 



accuracy [0]. We thus use this relation, instead of the charm quark mass itself, to test the 



theory. We note that Mf, and are calculated from the MS running masses, mb{Mb) and 
mc{Mc) using two loop QCD threshold corrections. In fact there are infra-red renormalons 
which make this perturbative evaluation ambiguous. However, these renormalon contribu- 
tions cancel in the mass difference. Moreover, it has been argued that this ambiguity cancels 
when extracting Vcb via semileptonic B decays p!7|][p!8| .p| 



Similarly, among the three light quarks there is one good relation, (reparametrization 
invariant and free of 0{mq) corrections) which severely constrains any theory of fermion 
masses. This is the Kaplan-Manohar-Leutwyler ellipse given by 

or 

rris 

Q = (3) 

mj. 

a 



'^The correction to the bottom mass {Mf, — mb(Aff,))/mf,(Aff,) at two loops is w 25% less than the re- 
sults found recently by Ball, Beneke and Braun using an infinite order resummation within a naive non- 
abelianization of the mass correction[^. The theoretical error for this result is also estimated to be about 
25%. By using two loop pole masses, combined with the mass relation for Mf, — we obtain a small value 



for mc{Mc). This value increases when using the scheme of BBB 18 
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where Q is the eUipse parameter. The experimental value for Q is obtained from a weighted 
average of lattice results and a chiral Lagrangian analysis, with significant contributions from 



the violation of Dashen's theorem ||19||. We follow Donoghue||20|| and assign a conservative 
10% error to the experimental value of . Note that md/rris derived from chiral Lagrangian 
analysis is not free of first order quark mass corrections, and hence a is much smaller for 
than for m^/ms |2^. In addition, we do not constrain mu/md independently, as suggested 
by Leutwyler [pT| , which requires additional input from an expansion in 1/Ncoior- Instead we 
quote values for mu/md as output. Finally we use given by the PDG[^. At the end of 
section V we consider the lower values for suggested by recent lattice calculations P5[]. 

The remaining parameters are more or less self evident. We just remark that the central 
value for Vd,-, as well as the error bars, has steadily decreased in the last 5 years, making it a 
very significant constraint .0 In addition, the value for Kib/Kfe changed dramatically in 1992. 
It changed from approximately 0.15 ± 0.05 to its present value 0.08 ± 0.02, where the errors 
were and continue to be dominated by theoretical model dependence. Clearly the systematic 
uncertainties were large but are now hopefully under control. 



3 Low energy observables — computed values 

In our analysis we consider the minimal supersymmetric standard model defined at a GUT 
scale Mg with, in all cases (but one), tree level GUT boundary conditions on gauge couplings 
and Yukawa matrices. In this one case, we include an arbitrary parameter, es, given by 

_ a^iMc) - ag ~ _ ^ ^ 

es = with ac = ai[MG) = a2[MG) (4) 

which parametrizes the one loop threshold correction to gauge coupling unification.]]] We 
also include 7 soft SUSY breaking parameters — an overall scalar mass mo for squarks and 
sleptons, a common gaugino mass M1/2, the parameters Aq, B and /x, and in addition we 
have allowed for non-universal Higgs masses, and m^/^.^ Thus the number of arbitrary 

^Note, the actual value of Vcb is correlated with the values of mb and mc jl^. We have not included this 
correlation in our analysis. 

^63 is calculable in any complete SUSY GUT. It is also constrained somewhat by the bounds on the 
nucleon lifetime 

^Note that if the messenger scale of SUSY breaking is Mpianck then our analysis is not completely self- 
consistent. In any complete SUSY GUT defined up to an effective cut-off scale M > Mq, the interactions 
above Mq will renormalize the soft breaking parameters. This will, in general, split the degeneracy of squark 
and slepton masses at Mq even if they are degenerate at M. On the other hand, bounds on flavor changing 
neutral current processes, severely constrain the magnitude of possible splitting. Thus these corrections must 
be small. In addition, in theories where SUSY breaking is mediated by gauge exchanges with a messenger 
scale below (but near) Mq, the present analysis is expected to apply unchanged, since in this case squarks 
and sleptons will be nearly degenerate at the messenger scale. The Higgs mases, on the other hand, are 
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parameters in the effective GUT includes the 3 gauge parameters, 7 soft SUSY breaking 
parameters and Uy Yukawa parameters. The number Uy and the form of the Yukawa matrices 
are model dependent. The models are discussed in sections IV and V. 

The effective theory between Mq and Mz is the MSSM. We use two loop SUSY renormal- 
ization group equations [RGE] for dimensionless parameters and one loop ROE for dimen- 
sionful parameters from Mq to Mz- However we have checked that the corrections to our 
results obtained by using two loop RGE for dimensionful parameters p6|| are insignificant. 



When in the framework of the MSSM we use Dimensional Reduction regularization in the 
Modified Minimal Subtraction Scheme. This renormalization procedure will be abbreviated 
as DR in the following text. 

When crossing the Mz scale we match the MSSM to the non-supersymmetric SU{?))c x 
f/(l)em gauge theory, i.e. our effective theory below the Z threshold is standard QCD and 
electromagnetism. In this effective theory we perform the calculations using Dimensional 
Regularization in the Modified Minimal Subtraction scheme. Quantities renormalized ac- 
cording to this prescription will appear with a standard subscript MS in the following 
formulas. Below the Z scale, three loop QCD renormalization group equations run quark 
masses and as down to lower scales. When this running crosses the h and c mass thresholds 
the number of flavors is reduced by one each time. Finally we end up with the three flavor 
QCD RGE as we arrive at the scale of 1 GeV , where the three light quark masses are eval- 
uated. In addition, QED running is included to one loop precision. (For further details see 



At Mz we calculate complete one loop MSSM corrections to the W and Z masses. Thus 
the W and Z pole masses are given by the formulae: 

Ml = ^glv' + 5Ml, (5) 
Ml = \{l9! + g!)v^ + 5Ml (6) 



where 6Ml and 5M| are the complete one loop self-energies |14l,[28| in the DR scheme and 



Qi, g2 are the DR gauge couplings evaluated at Mz in the MSSM. 

The Higgs vacuum expectation value v is an implicit function of soft SUSY breaking 
parameters and gauge and Yukawa couplings. It is determined by self consistently demanding 
minimization of the tree level Higgs potential. The actual value for v is found by minimizing 



The theoretical value of is calculated according to the formula 



probably dominated by new interactions which also generate a n term. It is thus plausible to expect the 
Higgs masses to be split and independent of squark and slepton masses. The parameter Aq could also be 
universal at the messenger scale. 
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The right-hand side of this equation is consistently evaluated in the DR scheme within the 
context of the MSSM. The electromagnetic structure constant a and sin^ Oy/ are determined 
by the relations among the gauge coupling constants ai and 



sin^ Ow{Mz) 



lai{Mz) + a2{Mz) 



a{Mz) = a2{Mz) sin'Ow^Mz) . (9) 

^ ^ contains the corrections from the W self-energy and vertex and box diagrams. 
We follow refs. (|T^,|^) closely, but not to all details: we don't include the vertex and box 



diagrams containing SUSY particles. These have a minor effect (see results in |T^) and 
we compensate for them with an additional 1/2% error on G^, as already mentioned in a 



previous footnote. In fact, we use a reverse procedure to the one introduced in [Q, where 
the equation (J^) is utilized to derive the MSSM value of sin^ 6w in the bottom-up approach 
starting from the precise value of G^. 

When crossing the Z threshold, gauge couplings are subject to two different kinds of 
corrections. That is to account for the states which are integrated out at this scale , and for 
the transition from the DR scheme to the MS scheme. First, when the SUSY particles, the 
top quark and the Z and W gauge bosons are integrated out , the threshold corrections to 
^^[dr ^^"^ ^[dr read 

1 ^ mg- Mg I Mt 

as 'n 12 ^ Mz Mz 3 Mz 

( all 
squarks) 

— = — ( - y bg ^ + - y bg ^ + - y og — ^ 

a 2tt ^ 9 tt ^ Mz 9 tt Mz 3 U ^ Mz 

( u — type { d — type ( charged 

squarks) squarks) sleptons) 

42 Mc, I Mh+ 16, Mt Mw, 

+ - > log H — log log 7 log 1 . (11) 

^ Mz 3 Mz 9 ^ Mz Mz' ^ ' 



The DR values of a^, (z = 1, 2, 3) are converted to the corresponding MS values using the 
relations 

, \, = , \, +C,; with Q = ^^^^, (12) 

where the quadratic Casimir C2{Gi)=N for SU(N) or equals zero for U(l). 

The strong coupling constant is now ready for comparison with the experimental value 
quoted in table [l|. The electromagnetic fine structure constant may stay within the DR 
scheme and it is corrected by an additional factor 0.0684 (|T^ and the references therein) 
to account for the running down to the zero momentum transfer value as extracted from 
Thomson scattering. 

Note that sin^^vF in eq.(^ represents the MSSM value, in the DR scheme, which can 



differ by a few per cent from the SM value - as found in [|IJ] . For completeness, we also 
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calculate the SM value, in the MS scheme, and quote both quantities in table |[ Our SM 
value of sin^ 6w{Mz) is obtained from the SM values of ai and a2 (by a relation analogous 
to eq.(D). The latter are derived from the corresponding MSSM values by applying the 
threshold corrections similar to ([T0| ) and (|ll|), with the W , Z and top quark terms absent. 



since these particles are kept in the SM. These are the leading logarithmic thresholds (LLT) 
which make up for the exclusion in the SM of the SUSY partners and non-minimal Higgs 
sector.0 A constant due to the change from the DR scheme to the MS scheme is added to 
5^2, according to eq.(|T2|). 

In the absence of direct observations, the SUSY sector of the MSSM is tested by all 
observables calculated at the loop level. Obviously, the strongest constraints come from the 
phenomena which originate directly in loop effects. Because of this reason pnew, as well as 
the process h 57, were added to the list of observables. Breaking of custodial SU{2) by 
physics beyond the SM is tested using 

new new\^J /i o\ 

M^ M| ' ^ ' 

where nJ^^(O) and nf^(O) stand for the contributions of physics beyond the SM to the 
vector boson self-energies at zero momentum ||16| . 

For the fermionic sector, we have at tree level m„ = Yu sin/3 nid = cos (3 and 
me = 1^6 cos /5 ^ , where rrij and Yi {% = u, d, e) are 3x3 mass and Yukawa matrices at Mz- 
In addition, at Mz we include the leading ( 0(tan/3)) one loop threshold corrections to the 
mass matrices [Q. The one- loop corrected fermion mass matrices are then diagonalized. The 
pole mass of the top quark is calculated from the diagonal running mass with the leading 
two loop QCD corrections included. We utilize the infrared fixed point behaviour of the 
top Yukawa coupling and the fact that as also changes only slightly between M^p and Mz, 
and evaluate the top pole mass from the relevant couplings at the Z scale.|3 The other 
quark and lepton masses are run down to their corresponding mass scales as explained in 
the beginning of this section. Mf, and Mc are evaluated using the two loop QCD corrections 
while the lepton pole masses are obtained with the help of the one loop QED corrections. 
Upon diagonalizing the one loop corrected mass matrices at Mz, we calculate the complete 

^Explicit forms of the LLT corrections to ai and a2 are given e.g. by the equations (17) and (18) in ref . [p^ . 
We agree with the authors of that the LLT approximation cannot proceed in a consistent way at this 
point since the SU{2) x U{1) gauge symmetry is broken at a scale comparable to the masses of sparticles 
being integrated out. The difference compared to the full calculation, however, is not very significant, and 
for our purposes will be neglected. Thus our SM value of sin^ 6\y contains a theoretical uncertainty of about 
1%, (which is larger than the experimental uncertainty) similar to the theoretical uncertainty assigned to 

^"^A more precise evaluation would use as{mt) and mt{mt) to obtain the pole mass Mt- The error obtained 
with our approximation, using values at Mz instead, partially cancels the one due to the non inclusion of 
the logarithmic corrections to the top mass, leading to a value which differs by less than 2 GeV from the 
actual value of Mf. 
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form of the CKM mixing matrix (note that the CKM matrix elements do not run below the 
Z scale), and then compute the theoretical value for needed to fit ei^^O|. The effects 
of SUSY box diagrams are neglected in this last calculation, while a potentially significant 
(0(tan/5)) one loop threshold correction to J |Q is included (automatically, since we start 
from the one loop corrected mass matrix 111^) as a leading higher order effect to the SM box 
diagram. The SUSY box diagram corrections can be as large as 15% and may affect 
the predictions for CP violating parameters, such as 77 or sin 2a and sin 2/3. We discuss the 
sensitivity of our results to these corrections in section V. 

Finally, special attention is paid to the branching ratio for h 57. We calculate the 
partial amplitudes for this process, as well as the amplitude for h s gluon (required to 
resum leading logarithmic corrections), using the complete quark and squark mixing matrices 



at low energies |3^. Large QCD corrections are resummed to leading order following 



The low scale in this evaluation is set equal to the pole mass Mb of the b quark, and the 
dependence of the branching ratio on this choice is neglected compared to the experimental 
error. 

We then form a function including the 20 low energy observables — 6 in the gauge 
sector {Mw, Mz, asM, Gfj,, as, p), 13 in the fermion mass sector (9 charged fermion masses, 
3 quark mixing angles and Bk) and the branching ratio for b ^ s'y. This function is 
then minimized self-consistently using Minuit by iteratively varying the GUT parameters 
with mo, Ml/2, and fi fixed. The procedure involves two nested minimizations. The GUT 
scale parameters mfj^ and ttlh^ are varied separately in the second nested minimization, 
for each step in the main minimization process where the rest of the free parameters are 
varied. Such a nested minimization is time consuming, but on the other hand it guarantees 
a high precision in the optimization process which is very sensitive to radiative electroweak 
symmetry breaking. Analytic tree level conditions for the potential minimum are used to 
determine v and tan/3. With these technical details we conclude the discussion on the 
theoretical calculation of the twenty observables given in table ^ 

In addition to the 20 observables discussed above, we have to deal with the experimental 
lower bounds on sparticle masses. We set lower limits of 30 GeV on squark and slepton 
masses, expecting significant positive higher order corrections to the squared mass which 
would raise it above the experimental lower bound of 45 GeV. We require that the light- 
est neutralino (if below the LEPl threshold) contribute to the invisible Z partial width as 
measured at LEPl - r{Z invisibles) = (—1.5 ± 2.7)MeV p5|]. We also demand that the 
lightest chargino be heavier than 65 GeV. Finally we require the pseudo-scalar Higgs, A, 
to have mass greater than 80 GeV. Note that these bounds are on the tree level running 
masses. In table |^ which shows predictions of the models discussed in the next sections, 
Higgs masses are displayed with the leading one loop corrections included.^ As can be seen, 

^^This CP odd state A, if light, (along with the other light Higgs states) receives significant corrections 
to its tree level mass in the regime of large tan /3 in which we work. If the pseudoscalar mass is less than 
Mz, it receives negative mass squared corrections; the same is true for the lightest CP even Higgs /i° [p6l. 
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the 80GeV lower limit on the tree masses provides enough room for the one loop corrected 
Higgs masses to stay above the experimental limits. In order to take all these bounds into 
account we have added a large penalty to the function whenever one of these bounds is 
exceeded, thus forcing the optimization procedure to stay within the experimentally allowed 
region in parameter space. 



4 Model 4 of ADHRS 

In this paper we have analyzed several models of fermion masses. We have studied model 4 
of ADHRS [0. In this model ny = 5. The model is defined by the following 4 operators in 
the effective theory at Mq- 



O33 = I63 lOi I63 (14) 

O23 = I62 lOi I63 
A A 

There are six possible choices for the 22 operator; all give the same 0:1:3 Clebsch 
relation between up quarks, down quarks and charged leptons responsible for the Georgi- 
Jarlskog relation | 



22 — 

(a) I62 4- 10i4^162 (15) 

OM A 

(b) I62 & lOi 4^ I62 
^ ^ A ' Sm 

A Ai 

(c) I62 — lOi I62 

<~>M <~>M 

(d) I62 lOi 4^ I62 

A 

(e) I62 lOi 4^ I62 

if) I62 lOi ^ I62 

Ai, A2 and A are adjoint scalars with vacuum expectation values [vev] in the B-L, Y or 
X (SU(5) invariant) directions of SO(IO) and Sm is a singlet with vev of order M. 
The resulting Yukawa matrices at Mq are given by - 
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where the Yukawa matrices are defined by the couphngs — HuQYuU + HaQY^D + HaLYeE 
and, for example, the fields Q{U) are left-handed Weyl spinors transforming as a weak 
doublet (singlet). 

The best fits give ~ 13 — 14 for 5 degrees of freedom. See figures la - c for contour 
lines of constant the mo — M1/2 plane for different values of Note that the regions 
with lowest contain significant (of order (4 - 9)%) one loop SUSY threshold corrections to 
fermion masses and mixing angles. These one loop SUSY threshold corrections scale roughly 
as ^^y^ or i^^. Hence, from figures la - c, we see that these corrections are necessary to 

mp rrig 

improve the agreement of the model with the data. 

In table ^ we give the results for the point labelled "II" in fig. la. All low energy 
observables are fit to within their 2a bounds. However the fits for 4 observables, [Vet, Vub/Vcb, 
Bk and Q) lie close to the 2cr border. Note that Minuit typically tries to equalize the 
contribution of all the observables to x^- We now argue that these results are true predictions 
of the theory. 

Consider the first three parameters. It was shown by Hall and RasinBH] that the relation 



(16) 



holds for any fermion mass texture in which the 11, 13 and 31 elements of the mass matrices 
are zero and perturbative diagonalization is permitted. Note X^, Ac are the up and charm 
quark Yukawa couplings evaluated at a common renormalization scale. A typical value for 
the right-hand side of the equation is 0.05 which is too small for the left-hand side by more 
than 20%. 

We now show that the fits for Vub/Vcb, Vcb and B^ are correlated. Consider the formula 
for ex given by 

ex ~ ( 77- VcbSin{^) ) Bk x (one loop factors) (17) 
V Vcb ) 

where the first factor is just the Jarlskog parameter, J - see eqn.(|]). We see that if VubjVcb 
is small, then Vcb and Bx must be increased to compensate. As a consequence, Vcb and Bx 

and mo arc GUT scale values, while /i is given at Mz- 
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are both too large. The addition of 13 and 31 mass terms modifies relation (|T6D and can, in 
principle, accomodate larger values of Vub/Vcb and thus lower values for Vcb and B^- 

Now consider the ellipse parameter Q. This parameter is strongly controlled by the 
Georgi-Jarlskog relation [p7| 

m, 1 m,, , , 

— ^ ^ 18 

vfid 9 me 

which is satisfied by model 4. This is an important zeroth order relation to try to satisfy. 
However unless there are small calculable corrections to this relation, it leads to values of 
ms/md ~ 25 and thus values of Q which are too large. Note, that introducing 13 and 31 
terms in the down quark and charged lepton mass matrices can also, in principle, perturb 
the zeroth order Georgi-Jarlskog relation. 

Thus the disagreement between model 4 and the data seems significant. It is unlikely 
that it can be fixed with the inclusion of small threshold corrections to the Yukawa relations 
at Before one adds new operators, however, it is worthwhile to consider the possibility 

that perhaps one of the experimental measurements is wrong. It is interesting to ask whether 
the agreement with the data can be significantly improved by removing one contribution to 
the function; essentially discarding one piece of dataQ In order to check this possibility 
we have artificially infiated the value of a in the analysis for several observables (one at a 
time), to see if this significantly improves the fit. Our results are given in figure 2, where we 
also state, in the figure caption, the preferred value of the observable with the infiated error. 
As one can see, we find no significant improvements with this procedure. 

On the other hand, we have found that we can indeed improve the results by adding 
one operator contributing to the 13 and 31 elements of the Yukawa matrices. We discuss 
this possibility in the next section. The additional terms correspond to one new effective 
mass operator. Of course there are many possible 13 operators. In this work we have not 
performed a search over all possible 13 operators. Instead we study two 13 operators which 
are motivated by two complete SO(IO) extensions of model 4. 

Finally, we should note that we have also studied the other 8 models of ADHRS. The 
values for these models are significantly larger than those for model 4 discussed here. 



5 Models 4(a,b,c) of LR 

In this section we analyze two models derived from complete SO(IO) SUSY GUTs discussed 
recently by Lucas and S.R.[0. The models were constructed as simple extensions of model 4 
of ADHRS. The label (a,b,c) refers to the different possible 22 operators which give identical 

^•^Note, there are 2 types of threshold corrections. First, there are the higher dimension operator corrections 
to the 4 operators; obtained when integrating out the Froggatt-Nielsen massive intermediates. These are 
expected to be of order 10%. Second, there are the usual higher loop corrections. 

^''We thank M. Barnett for bringing this idea to our attention. 
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Clebsch relations for the 22 element of the Yukawa matrices. However in the extension to 
a complete SUSY GUT these different operators lead to inequivalent theories. The different 
theories are defined by the inequivalent U(l) quantum numbers of the states. When one 
demands "naturalness", i.e. includes all terms in the superspace potential consistent with 
the symmetries of the theory one finds an additional 13 operator for models 4a and 4c given 
by- 



O 



13 



(19) 




(a) 
(c) 



Model 4b, on the other hand, is identical to model 4 of ADHRS when states with mass 
greater than Mq are integrated out. The results for this model are identical to those pre- 
sented in the previous section; thus we will not discuss it further. The other two models 
include one new effective operator at the GUT scale. The addition of this 13 operator intro- 
duces two new real parameters in the Yukawa matrices at Mq] thus we have ny = 7. Models 
4a and 4c differ only by the 13 operator. The resulting Yukawa matrices are given by - 



where the new Clebsches for models 4(a,c) are given in table ^ 

We find that model 4a gives a best fit ~ 4 for 3 dof, while model 4c gives ^ 1 
for 3 dof. For model 4c, our best fit with = 0.168 is found for /i = IGOGeV, M1/2 = 
400GeV,mo = 2900GeV^. Our results are presented in figures [3(a - c)] for model 4a and 
figures [4(a - c)] for model 4c. 

In order to understand why the 13 operator resolves the problems discussed in the pre- 
vious section, consider the approximate formulas (valid to 10% and particular to model 4c) 

""^^Note, models d, e and f have the second family I62 coupled directly to lOi and a heavy 16. If this 
coupling is as large as the third generation Yukawa coupling, then we would obtain excessively large flavor 
changing neutral current processes, such as /i — > ej. Thus these models were not considered in [Q. 
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for Vub/Vcb, Vcb, As, Ad, A^, Ag where the latter are the diagonahzed Yukawa eigenvalues and 
all are evaluated at the GUT scale. We have 

Vub/V^b = ^ (20) 



AB 



cb 



9A 
Xs = E 

729C^f 31BD s 

A^ = 3E 

243 C / 109 
K = — — 1 - ■ e 



E V 27 C A 



(21) 



Recall, the contribution of the 13 operator is proportional to D. Consider first the electron 
and down quark Yukawa couplings. The electron gets a significant correction from the 13 
operator whereas the correction to the down quark is much smaller. When one adds the 
13 operator to a model which is already fit to the lepton masses, one must readjust the 
parameters. In order to keep the electron mass fixed, one needs to increase the parameter 
C [note, in the fits 6 ~ 27r]; the parameter E is kept fixed in order not to change the muon 
mass. As a consequence the ratio Vub/Vcb increases, which allows for smaller values of Bk- 
In addition we see that the Georgi - Jarlskog relation is now given by 

1 _ 109 BD_ -i 6 \ _ 
^ 27 CA*^ 1 




As a consequence, the ratio 1/9 is further decreased; hence Q is reduced. Thus one new 
operator is able to resolve four problems. 

Let us now discuss our results for model 4c. We find the preferred region of SUSY 
parameter space with a fixed value of /i = 80GeV corresponds to M1/2 > 220GeV and 
rriQ > 300GeV. The lower bound on mo depends slowly on M1/2 while the lower bound on 
Ml/2 appears to be independent of mo. These bounds are rather distinct (the value rises 
steeply when getting closer to these values) since they result from the lower experimental 
limits on sparticle masses. No observable from our list in table |1| places a decisive constraint 
which would exclude a portion of the (mo, M1/2) parameter space at fixed fi = 80GeV. In 
other words, for this value of fi everything in the (mo, M1/2) plane that has not been excluded 
by direct experimental searches, is allowed (provided the other parameters are subject to the 
optimization procedure). 

As /i increases (see figs. 4b and 4c for ^ = 160, 240 GeV, resp.) the lower bound 
on Ml/2 goes down since it is determined by the LEP limits on the masses of the lightest 
chargino and neutralino. The latter become proportional to, but less than, M1/2 in the limit 
of large fi. In addition, as /i increases the profile in the mo direction starts to change 
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more smoothly, and the contour hnes of constant shift more towards the higher values 
of tuq. This is explained by the fact that the constraints placed by the observables which 
receive one loop corrections enhanced by tan/3 start to be significant. (Especially the b quark 
mass correction and subsequently the correction to Mb — M^.) Since these one loop SUSY 
threshold corrections scale roughly as yuMi/2/mQ or /iAss/mQ we see that with increasing 
a larger suppression by uiq is required to keep the potentially large corrections under control. 
The effect clearly shows up in figures 4a,b,c combined. 

Our results show that the SUSY sector can make its presence visible in the analysis of 
the SM fermion mass parameters, especially those suppressed by tan/3 at tree level, through 
the low energy threshold corrections. For ~ 1 the latter are up to 9%. Two notes are 
in order at this point: First, note that these corrections have the potential to be as large 
as 30-70% 00. It would be interesting if a window in parameter space with such large 
SUSY contributions was favored. However, our numerical optimization never disclosed such 
a region. Second, it is interesting to note that the best fit is found in the region of very 
large mo, where the effect of the SUSY threshold corrections to fermion masses and mixings 
is, in fact, minimized. As a result, in this region, with < 1 and the SUSY corrections 
to fermion masses negligible, the effective number of degrees of freedom is actually larger 
than 3, since there are 7 parameters in the Yukawa matrices determining the 13 low energy 
observables in the fermion mass sector. This means that the Yukawa sector of the selected 
model does actually a much better job than appears at first glance. 

For the opposite sign of /i there is a significant restriction from the process b — > 57. As, 
has been noticed by many authors, the SM amplitude with the W in the loop is sufficiently 
big to explain the measured branching ratio. In the large tan/3 regime of the MSSM however, 
there are substantial partial amplitudes due to the charged Higgs and chargino loops. The 
charged Higgs amplitude is always of the same sign as the SM amplitude, while the sign 
of the chargino amplitude depends on the sign of ixAt. In order to get an agreement with 
the measured rate, the chargino amplitude, especially with large tan/3, has to enter with the 
opposite sign, and as a consequence a specific sign of ^At is preferred. Due to the boundary 
conditions at the GUT scale which are used throughout this paper we always run into the 
region at the Z scale where At < 0, and as a result the sign of /i is fixed. (;U > 0, in our 
convention.) With the bad sign of /i, the best values are about 30.|^ For this reason 
we do not quote results or figures for /i < 0. Interestingly enough, even if we exclude the 
branching ratio for b s'~f from the list of observables and repeat the optimization procedure 
for n < 0, the conspiracy between the one loop corrections proportional to /i (see the two 

Clearly, such large corrections would not violate perturbativity, since the main reason for them being 
large would be that the corresponding tree level values happen to be tan/3 suppressed. That means that the 
higher order corrections would be well under control as for any other quantity not suppressed at tree level. 

^^Obviously, if superpartners are very heavy, the SUSY spectrum decouples and the above statement is 
incorrect. However, for large tan/?, ^^very heavy" in this case means deep in the TeV region which we don't 
consider in our analysis. 
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previous paragraphs) does not work very well and significantly worse values are obtained 
than in the case /i > 0. 

Finally note that less than about 65GeV is excluded since there would inevitably be a 
light higgsino-like neutralino and chargino; already ruled out by experiment. Moreover, the 
large values of mo selected by our fit, imply a suppression of the supersymmetric corrections 
to the anomalous magnetic moment of the muon, which could otherwise lead to very strong 



constraints on the allowed parameter space in the large tan /3 region 



In tables ^ - |^ we give the values of the initial parameters which minimize for fixed 
particular values of /x, mo, M1/2 labelled as points I-III marked in figures 4a-c. For each 
low energy observable we quote its computed value, the partial contribution to and the 
relative magnitude of the SUSY threshold corrections [in %]. We also include in tables ^ 
and H values for the SUSY and Higgs spectra, the unitarity triangle parameters p, r] together 
with the corresponding angles a, /3 and 7, measurable in neutral B decays, and the predicted 
ratio for m^/md. 

Note that the values for a, /3 and 7 do not significantly change across different points in 
the SUSY parameter space. This fact is explicitly evident from figure 5. Each "x" -symbol 
in this figure stands for a point anywhere in the SUSY parameter space with the minimum 
of less than 3, with 3 degrees of freedom present both in model 4c and model 4a. Note 
that, under such a restriction in model 4a, only the corner in the SUSY parameter space 
with very large mo and M1/2 contributes to figure 5b. Nevertheless, these two models give 
similar, narrowly spread predictions for sin2Q; and sin2/? , grouped around the values 0.95 
and 0.52, respectively, while a much larger region is allowed for these two observables as a 



result of a general Standard Model analysis 44 



We have also checked the sensitivity of our results, for point I in fig. 4a, to a potential 
10% positive enhancement of e^- We have found that the total stays the same to within 
a few per cent and that our predictions for sin2a (sin2/3) change from 0.953 (0.513) to 0.957 
(0.527). Indeed, over a large range of points, we checked that the predictions for sin2a 
(sin2/5) change by less than 1% (3%). Thus our results, for model 4c, are not very sensitive 
to SUSY box corrections to ex- On the other hand, for point II (fig. 3a) in model 4a, 
we find that a 10% positive enhancement of significantly improves the agreement with 
experiment; x^ changes from 6.99 to 4.57. This is because in this case, -Ba' was a significant 
constraint on the model. 

It is an interesting feature of most of the points in the SUSY parameter space that the 
best fits favor a rather large negative GUT threshold correction €3 (eqn. In figure 6a the 
sensitivity of the best fits to different lower bounds on €3 is displayed for values of \i and 
Ml/2 fixed and mo varied. In addition table |iy shows how the selected GUT scale parameters 
Mq and A (the 33 element in the Yukawa matrices) vary in the best fits at mo=700GeV 
(point II), in order to compensate for changes in the lower bound on 63. Our study shows 
that even small positive GUT threshold corrections 1-2% to are plausible in model 4c. 
As one can see in table 10, this would require a rather low value of the GUT scale, below 
1 X lO^^GeV, and lower ac- As a result, the best fit value of as{Mz) ~ 0.116 is getting closer 
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to its central value 0.118, while is getting worse. On the other hand, the A parameter of 
the Yukawa matrices is quite insensitive to the change in and to the induced changes in 
ac and Mq, and so are the masses of the heaviest generation. 

We note that our results in this regard differ from those of a recent paper by Pierce 



et al.[0 in which a general MSSM analysis is performed. We believe that a significant 
part of the difference results from our assumed theoretical uncertainties in the electroweak 
observables. Recall, in our analysis we have assumed an overall 0.5% theoretical uncertainty 
(a) for Mz, Mw, and a em and a 1% uncertainty for G^. As commented previously, we believe 
these uncertainties correctly account for the theoretical errors introduced by neglecting higher 
orders in perturbation theory and by the numerical analysis. However for comparison, we 
have changed the theoretical uncertainties for Mz, (Xem, and G^ to 0.1% and for to 
the experimental value of 130 MeV; the new results, for model 4c, are displayed in figure 6b 



and table |TT|. In table |Tl| we see that in this case as varies significantly with changes in e^. 
Moreover, as increases the quality of the fit quickly deteriorates. This is seen graphically 
in figure 6b. We find, for point II in model 4c with eg = -3% (-2%), = 7.49 (13.97) with 
the dominant contributions to coming from Mz - 2.11 (3.67), - 0.89 (1.81), uem ^ 
0.74 (1.46) and nib^Mb) - 2.05 (3.69). Thus the theoretical uncertainties for the precisely 
measured electroweak parameters allow for quite a bit of flexibility in the values of eg and 
asiMz). 

Since we have focused on fermion masses we have not discussed the Higgs sector of 
particular models in great detail in this work. However, in figure 7 we show the sensitivity 
of the best fits to the assumed lower bound of the pseudoscalar mass, at fixed values of 
/i=80GeV and Mi/2=240GeV in model 4c. As one can see, requiring a heavier pseudoscalar 
(A) forces the model into the region of larger tuq and thus heavier squarks and sleptons. Or 
in other words, in order to have light squarks and sleptons, one necessarily also has a light A. 

To understand these properties, it is instructive to analyze the renormalization group 
evolution of the mass parameters. For a top quark mass Mt ^175 GeV, the low energy soft 
supersymmetry breaking mass parameters are approximately given by p3 



^ ^(m^JO)-m2)-2.5Mi%, 



m 



\ (mi^(O) - mg) - 2.3Mi% - ^ (2 + m|^(0)) , 



, 5 mfr (0) + m?r fO) 



„ 5 2 m?r fO) „ , , 

ml ^ |il^ + 5M,%, (22) 

where mg and m/j are the left and right handed sbottom mass parameters, and we have ig- 
nored the Aq contribution, as well as the small rnnS^iid) effect on the running of muSvau^. 
Moreover, in the large tan/? regime there is a simple relation between the Higgs mass pa- 
rameters and the pseudoscalar Higgs mass, namely [p^ , 

m\ ~ m\^ - - Ml. (23) 
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Using the fact that mjj^ ~ — M|/2 — /i^, we obtain, 

m^jO) - ^ 4M^, - ^ (/.^ + ^) . (24) 

Since the best fit is obtained for moderate values of the one-loop down quark mass correc- 
tions, low values of the fi parameter /i^ -C ttiq are preferred. For M1/2 > ~ 0{Mz), an 
approximate relation between m//„(0) and mo is obtained as a function of the gaugino mass 
parameter M1/2. Moreover, 

0.9 mi^(O) - mi^(O) ^ 0.1 (2 - 3 M^^) + ^(m^ + M|). (25) 

From the above equations, the qualitative behaviour of our solutions may be understood. 
For instance, for mo — 3Mi/2, as is the case at point II on fig. 4a (see table V), and 
Ml/2 > f^,Mz, from Eq. (|2^) we get that rriH^iO) ^ 1.2 tuq. Moreover from Eq. (^), we 
obtain that mi^^(O) may vary from values of order 1.4 mo, for low values of m^, up to values 
of order 2 mo for larger values of m^, without inducing a problem in the sbottom sector. 
Larger values of mj^^(O), however, induce lower values of the sbottom mass parameters, 
implying unwanted large values for the gluino-induced bottom mass corrections. Therefore, 
the best fit is obtained for values of the CP-odd Higgs mass close to its experimental bound, 
or, equivalently, for values of m/^^(0) of the order of 1.4 mo. Larger (smaller) values of 
mo/Mi/2 imply smaller (larger) values of m//^^(0)/mo, as is clearly seen from Eqs. (|2^), 
(^) (see Tables III- VII). 

The relation between the CP-odd mass and mo may also be understood from the above 
equations. For large values of the CP-odd mass, and mo > M1/2, ''^h^(O) must be significantly 
larger than fnjj^{0). For these large values of m'^^{0), large values of mg are required in order 
to avoid a very low value of the sbottom mass. The results displayed in fig. 7 are just a 
reflection of this fact. 

To conclude this section let's discuss how the performance of models 4a-c changes if recent 
lattice results on the low value of the strange quark mass prove to be correct With 
ms{lGeV) = (120 ± 25)MeV replacing the value (180 ± 50)MeV of table |1|, the value 
of model 4c at point II changes from 0.731/3dof to a significantly worse value of 5.44/3dof. 
More than a half of it comes from the rUg value itself, 163.4MeV at the best fit, indicating 
that model 4c cannot get to such low values of the strange mass as reported by the lattice 
groups, even with the help of the fifth effective operator. The rest of the value is being 
shared by a low value for as{Mz), a high Vcb and Bk, and low Vub/Vcb, contributing by less 
than 0.8 to the total each. We checked that this behaviour of model 4c is typical for 
different points as well, and that the best fit values of the strange quark mass m^ < 160MeV 
never occurs. Model 4a, on the contrary, favors values of ms{lGeV) less than 180MeV. 
However, decreasing m^ cannot cure the problems of this model: high values of Bk and Vcb, 
and a low value for Vub/Vcb- At the same point II, these observables contribute 2.85, 0.75 
and 1.76, respectively, to which totals 6.82 despite values of m^ as low as 136.5MeV. This 
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represents only a marginal improvement to the value of 6.99 at the same point with the 
original larger experimental value of the s-quark mass. Finally, the performance of model 
4 gets slightly worse with the low value of m^: total changes from 14.01 to 16.16 as 
goes from 163MeV down to 156MeV in the best fits at point II. 



6 Discussion and Conclusions 

The results of our analysis, as well as the whole project presented here, may be understood 
from two different perspectives. 

The emphasis of this paper is to reanalyze the fermion mass sector in the context of the 
minimal SO(IO) SUSY GUT in order to gain insight into the underlying flavor physics. From 
this perspective, we analyzed the best working model of ADHRS with four effective operators 
in the Yukawa sector at the GUT scale. The best fit allows us to assign a confidence level 
of about 3-4% to this model. Next, we found that the addition of a new operator, giving 
rise to the 13 (and 31) entries in the Yukawa matrices, may improve the performance of 
the model. Substantial improvement however, is not automatic, as has been evidenced with 
model 4a that has a confidence level < 50% for the best fits in the corner of SUSY parameter 
space with very large mo and M1/2, and large /i. On the other hand, we showed that model 
4c provides an excellent fit to all 20 low energy observables, with confidence level better 
than 68% in a large region of the allowed SUSY parameter space. We have not performed 
a complete search over all possible 13 operators. Model 4c was picked as a candidate model 
suggested by a recent formulation of a complete SO(IO) SUSY GUT|^. Whether or not this 
particular model is close to the path nature has chosen remains to be seen. One important 
test will be via the CP violating decays of the B. We predict a value for sin2a which is 
insensitive to the SUSY breaking parameters (see figure 5 and tables ^ and P), whereas in 
the SM the value of sin2a is unrestricted||44||. Another important test will come from nucleon 



decay rates which is discussed in a recent paper . 

As a separate matter, we considered the sensitivity of our results to different lower bounds 
for €3. For the assumed theoretical uncertainties, as given in table ||, our analysis of model 4c 
indicates that the best fits favor a negative GUT threshold correction to of the order -(4- 
5)%. We have checked for the whole mg - M1/2 plane at fixed /i = 80Gel^ that a correction 
of -3% instead of -(4-5)% can easily be accomodated and gives basically the same best values 
of as obtained in the unrestricted analysis. In this case, for 63 = —3%, the contour lines 
of figure 4a move almost uniformly towards higher values of mo by about 50 GeV. Note that 
the study of the complete SO(IO) model0] giving rise to the effective operators of model 
4c shows that such negative corrections to ^^(Mg) can easily be obtained without any fine 
tuning. To understand how sensitive our results might be to this threshold correction we 
studied a narrower region in the parameter space, determined by fixed Mi/2=240GeV at the 



previously fixed /i=80GeV (see figure 6a and table [1^). It turns out that €3 > remains a 



possibility, at the expense of lower GUT scale and ac-, resulting in a worse agreement with 
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the low energy value of G^. Fermion masses are not very sensitive to these variations. In 
addition, from table |10| it is clear that as increases, so does ^^(M^). These results are 
however strongly dependent on the magnitude of our assumed theoretical uncertainties as 



discussed in the previous section and shown in figure 6b and table 11 



Note that model 4c is also flexible enough to accomodate recent analyses of the Rb 
anomally and low values for as{Mz) 0. For the Higgs sector of the MSSM, the best 
fits at a number of points (e.g. also at the sample points I - III; see tables ^ and |^) end up with 
low values of the pseudoscalar Higgs mass rriA- This preference is not strong though, as can 
be inferred from figure 7 where we show how the quality of the fits change with the increasing 
lower bound on m^. However, the low value of itia-, such as at point III2 for instance, leaves 
the door open for a natural explanation of the l-2a increase in the partial width Z bb, 
without ever asking for it in the course of our analysis. At the same time, as can be argued 
on general grounds, one of the CP even Higgs states is always almost degenerate with the 
pseudoscalar state. For niA < Mz it is the lighter state, while for > Mz it is the heavy 
Higgs state H^. In the latter case, the mass equals M^cos^2/3 + rad. corrections where 
the radiative corrections are rather substantial in our case. Thus our study indicates that 
there are no additional constraints on a light CP even Higgs state between the experimental 
limit of about 55GeV (for tan f3 ~ 60) and the upper MSSM limit of about 130GeV. 

Another interesting feature of the best fits is the tendency towards lower values of 
as{Mz) ~ 0.115. As with low m^, it just happens to be in agreement with a potentially 
significant positive pseudoscalar contribution to Ri,. The origin of rather low in the best 
fits can be traced to the optimization efforts to suppress the gluino correction to an already 
large enough value for mb{Mz)- This correction has to be positive, since its sign is correlated 
with the sign of the chargino-stop partial amplitude to the process 6 — 57 , and the latter 
is fixed in the large tan/3 regime, as described in the previous section. At the same time 
lower values of as{Mz) are welcome in order to suppress the QCD logarithmic corrections to 
b S7. Here the point is that with tan/3 large the chargino-stop loop tends to be too large, 
most of the times outweighing by too much the W and loops which contribute with the 
opposite sign. Thus the optimization favors lower as{Mz) as the way to make up for a rather 
large net b ^ s'y amplitude by minimizing the enhancement from the renormalization to the 
Mb scale. 

Our analysis can also be viewed from a different perspective. Let us neglect for a moment 
the underlying GUT physics and the origin of the Yukawa matrices at the GUT scale and 
view this analysis simply as an MSSM global fit in the large tan/? regime. (See also ||48|| .) 
/^From this new perspective a special feature of our approach is that we run complete sets 
of fermion and sfermion mass matrices down to the Z scale, instead of just the leading 33 
elements. We now try to draw some general conclusions coming from such an analysis. 

Our results suggest that there is no narrow, strongly preferred region in the SUSY 
(mo, M1/2) parameter plane for low values of the /i parameter /i(M^) ~ 80GeV^ . Hence, at 
present one cannot make strong conclusions in the large tan/5 regime about the masses of 
squarks, sleptons and gluinos, which leaves open various channels for Tevatron and LEPII 
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experiments. 

In addition, a definite statement can be made about the sign of /i following from the 
structure of the partial amplitudes to the process b — > 57, as already explained in more 
detail in the previous section. In the conventions we use fi has to be positive. Note, this 
conclusion, however, depends on the assumption of universal squark and slepton masses and 
universal A parameter at Mq which constrains the sign of At at Mz- 

Finally, by keeping the complete 3x3 mass matrices for both fermions and sfermions 
one can study flavor dependent processes in a theory which fits the low energy data; for 
example, in rare B and K decays, B — B mixing, lepton flavor violating processes or even 
Z bb. 

In summary, we performed a detailed global analysis of independent SM parameters 
(plus the extra low energy data constraining the SUSY sector) in several models based on 
SO (10) SUSY GUTs. At the present time, when direct evidence for physics beyond the 
Standard Model evades experimental observations, this kind of analysis serves as the best 
test of new physics and actually starts to compete with SM electroweak precision tests. 
Our ultimate goal is to identify a set of effective theories, defined at a scale M (M = 
M Planck, M string OT Mqut) which accuratcly fit the low energy data and thus determine 
the boundary conditions at the scale M for some more fundamental theory valid above M. 
We have shown that this set is non-trivial; model 4c is one element of the set. More elements 
of the set need to be found in order to determine (a) whether they can be distinguished by 
purely low energy measurements, or (b) whether one or more elements can be constructed 
as the low energy limit of a string. 

Finally, some theoretical uncertainties in the present analysis can and should be removed 
in the future. By including the SUSY box contributions to we can remove the additional 
1/2% uncertainty included in the evaluation of for this quantity. In addition, the SUSY 
box corrections to ex can be as large as 15%||3^. These should be included in order to remove 
the uncertainties in the predictions for CP violating parameters; although these uncertainties 
do not appear to be significant at the representative points studied in model 4c. 
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Table 1: Experimental observables 
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Central value 


a 


Mz 


91.186 


0.46 


Mw 


80.356 


0.40 




1.166 ■ 10"^ 


1.2 ■ 10-^ 




137.04 


0.69 


as{Mz) 


0.118 


0.005 


Pnew 


-0.6 • 10-3 


2.6 ■ 10-3 


Mt 


175.0 


6.0 


mi,{Mb) 


4.26 


0.11 


Mh-M, 


3.4 


0.2 


rus 


180 


50 




0.05 


0.015 




0.00203 


0.00020 


Mr 


1.777 


0.0089 


M, 


105.66 


0.53 


Me 


0.5110 


0.0026 


^us 


0.2205 


0.0026 


Vch 


0.0392 


0.003 


Vub/Vcb 


0.08 


0.02 


Bk 


0.8 


0.1 


B{b^ 57) 


2.32 ■ 10"^ 


0.92 ■ 10"^ 
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Table 2: u Clebsches for models 4(a) and 4(c) 



model 


Uu 


< 








< 


a 


-4/3 


1/3 


-2 


-9 


-54 


3 


c 


-4/3 


1/3 


2/3 


-9 


-54 


-1 



Table 3: Model 4 - Results at point II (see point II on fig. la) 
Initial parameters: l/ac = 24.33, Mg = 3.29-lOi^GeV, eg = -4.72% , 

A = 0.734, B = 5.67■10-^ C = 0.865■10-^ E = 1.04-10-2, $ = 1.76, 
H = SO.OOGeV, mo = 700.00GeV, M1/2 = 240.00GcV, 
mnj-mo = 1.39, niH^/mo = 1.22, Aq = 387.79GeV, =120.01GeV2 



Observable 


Computed value 


Contribution to 


SUSY corrections 


[in% 




Mz 


91.12 


< 0.5 




Mw 


80.37 


< 0.5 






1.166-10-^ 


< 0.5 




O^EM 


137.0 


< 0.5 


1.44 


Cis{Mz) 


0.1162 


< 0.5 


12.93 


Pnew 


+ 1.74- 10-^ 


< 0.5 




Mt 


173.9 


< 0.5 


0.75 


mb{Mh) 


4.360 


0.82 


5.14 




3.146 


1.61 


7.90 


rus 


162.6 


< 0.5 


3.68 


md/rris 


0.0461 


< 0.5 


0.00 




0.00173 


2.19 


1.66 


Mr 


1.777 


< 0.5 


-1.94 




105.6 


< 0.5 


-1.46 


Me 


0.5113 


< 0.5 


-1.46 


Vus 


0.2215 


< 0.5 


0.00 


Kb 


0.0450 


3.79 


1.53 


Kfe/Kfe 


0.0463 


2.84 


0.00 


Bk 


0.9150 


2.10 


-3.09 


B{h 57) 


2.388 • 10"^ 


< 0.5 




TOTAL 


14.012 
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Table 4: Model 4c - Results at point I (see point I on fig. 4a) 

Initial parameters: I/cig = 24.43, Mg = 2.50- lO^^GeV, eg = -4.76% , 

A = 0.764, B = 5.26-10-^ C = 1.10■10-^ D = 4.63■10-^ 6 = 5.70, E = 1.25-10-2, $ = 1. 

= SO.OOGeV, mo = 400.00GeV, My2 = 280.00GeV, 
majmo = 1.77, niH^/mo = 1.59, Aq = 322.21GeV, Bfx = 120.00GeV2 



W Uot-i VclUlLy 




V^UiililUUllUii lU ^ 


SUSY corrections 


[in%] 




IVlz 


Q1 19 








80 34 








1.164 • 10-^ 


< 0.1 






136.9 


< 0.1 


1.25 


as{Mz) 


0.1132 


0.93 


12.66 


Pnew 


+9.87-10-5 


< 0.1 




Mt 


173.5 


< 0.1 


0.93 




4.311 


0.22 


4.96 


Mh-M, 


3.499 


0.24 


6.64 


rus 


184.6 


< 0.1 


4.95 


md/rus 


0.0496 


< 0.1 


0.00 




0.00205 


< 0.1 


2.01 


Mr 


1.776 


< 0.1 


-4.00 


M, 


105.7 


< 0.1 


-2.76 


Me 


0.5110 


< 0.1 


-2.76 




0.2205 


< 0.1 


0.00 




0.0406 


0.21 


2.49 




0.0744 


< 0.1 


0.00 


Bk 


0.8196 


< 0.1 


-5.05 


B{h S7) 


2.575 - 10-^ 


< 0.1 




TOTAL 


2.039 
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Table 5: Model 4c - Results at point II (see point II on fig. 4a) 

Initial parameters: Ijac = 24.36, Mq = 3.17-lO^^GeV, eg = -4.89% , 

A = 0.807, B = 5.44-10-2, C = LIS-IO^^ D = 4.94■10-^ 6 = 5.71, E = l.Sl-lO-^, $ = 1. 

H = 80.00GeV, mo = 700.00GcV, M1/2 = 240.00GcV, 
majmo = 1.42, niH^/mo = 1.24, Aq = 458.35GeV, Bfx = 120.66GeV2 



W Uot-i VclUlLy 




r~^/~\TTt"nT^n "t" 1 \'r\ '\/^ 
V^UiililUUllUii lU ^ 


SUSY corrections 


[in%] 




IVlz 




09 






80 38 








1.166-10-^ 


< 0.02 






137.0 


< 0.02 


1.43 


as{Mz) 


0.1151 


0.34 


12.78 


Pnew 


+1.87-10-^ 


0.09 




Mt 


175.7 


< 0.02 


0.74 




4.287 


0.06 


5.43 


Mb-M, 


3.440 


0.04 


7.56 


rus 


189.0 


0.03 


3.68 


md/rus 


0.0502 


< 0.02 


0.00 




0.00204 


< 0.02 


1.78 


Mr 


1.776 


< 0.02 


-2.08 


M, 


105.7 


< 0.02 


-1.50 


Me 


0.5110 


< 0.02 


-1.50 




0.2205 


< 0.02 


0.00 




0.0400 


0.07 


1.58 




0.0772 


< 0.02 


0.00 


Bk 


0.8140 


< 0.02 


-3.18 


B{h 57) 


2.382 • 10-^ 


< 0.02 




TOTAL 


0.7306 
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Table 6: Model 4c - Results at point IIIl (see point IIIl on fig. 4b) 
Initial parameters: Ijac = 24.51, Mq = 3.33-lO^^GeV, eg = -4.34% , 
A = 0.852, B = 5.63-10-^ C = 1.21■10-^ D = 5.06■10-^ 5 = 5.70, E = 1.36-10-2, $ = 1.02, 

H = leO.OOGeV, mo = MOO.OOGeV, My2 = 170.00GeV, 
majmo = 1.33, niH^/mo = 1.14, Aq = -982.43GeV, Bfx = 123.72GeV2 



W Uot-i VclUlLy 




r~^/~\TTt"nT^n "t" 1 \'r\ '\/^ 
V^UiililUUllUii lU ^ 


SUSY corrections 


[in%] 




IVlz 




09 






80 38 








1.166-10-^ 


< 0.02 






137.0 


< 0.02 


1.79 


as{Mz) 


0.1154 


0.26 


12.89 


Pnew 


+ 1.50- 10-^ 


0.08 




Mt 


176.9 


0.10 


0.62 




4.275 


< 0.02 


6.27 


Mh-M, 


3.429 


0.02 


8.76 


rus 


188.3 


0.03 


2.85 


md/rus 


0.0512 


< 0.02 


0.00 




0.00203 


< 0.02 


1.68 


Mr 


1.777 


< 0.02 


-1.41 


M, 


105.7 


< 0.02 


-0.82 


Me 


0.5110 


< 0.02 


-0.81 




0.2205 


< 0.02 


0.01 




0.0397 


0.03 


1.62 




0.0797 


< 0.02 


0.00 


Bk 


0.8000 


< 0.02 


-3.26 


B{h S7) 


2.311 • 10-^ 


< 0.02 




TOTAL 


0.5868 
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Table 7: Model 4c - Results at point III2 (see point III2 on fig. 4c) 
Initial parameters: l/a^ = 24.65, Mq = 2.88-lO^^GeV, eg = -4.45% , 
A = 0.888, B = 5.89-10-2, C = 1.23■10-^ D = 5.69■10-^ 5 = 5.74, E = 1.40-10-2, $ = 1.02, 

H = 240.00GeV, mo = MOO.OOGeV, My2 = 170.00GeV, 
majmo = 1.33, niH^/mo = 1.14, Aq = -1079.39GeV, Bfx = 128.50GeV2 



W Uot-i VclUlLy 




V^UiililUUllUii lU ^ 


SUSY corrections 


[in%] 




IVlz 












<^ 10 






1.165 - 10-^ 


< 0.10 






136.9 


< 0.10 


1.85 


as{Mz) 


0.1124 


1.25 


12.53 


Pnew 


+4.60- 10-5 


< 0.10 




Mt 


176.5 


< 0.10 


0.61 




4.306 


0.18 


9.73 


Mh-M, 


3.486 


0.19 


13.45 


rus 


178.5 


< 0.10 


4.09 


md/rus 


0.0500 


< 0.10 


0.00 




0.00203 


< 0.10 


1.91 


Mr 


1.776 


< 0.10 


-1.93 


M, 


105.7 


< 0.10 


-1.06 


Me 


0.5110 


< 0.10 


-1.06 




0.2205 


< 0.10 


0.01 




0.0402 


0.11 


1.73 




0.0777 


< 0.10 


0.00 


Bk 


0.7949 


< 0.10 


-3.49 


B{h S7) 


2.269 - 10-^ 


< 0.10 




TOTAL 


1.9477 
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Table 8: Model 4c - More Results at points I and II (see points I and II on fig. 4a) 
In Higgs sector, we quote all masses with one loop corrections included. For superpartners, 
masses are at tree level. 

For squark and slepton masses, the first two columns are for the third family which are sig- 
nificantly split, while the third and fourth columns are mean values for the nearly degenerate 
states of the second and first families, respectively. 



Observable 


Predictions at point I 


Predictions at point II 


mu/md \iGeV 


0.409 


0.440 




0.2342 








0.2334 








0.2323 








0.2320 






sin 2a 


0.953 








0.966 






sin 2/3 


0.513 








0.522 






sin 7 


0.935 








0.928 






P 


-0.120 








-0.130 






V 


0.316 








0.325 






tan/3 


52.77 








54.38 






CP even Higgses 


73.43 


110.81 






75.24 


110.35 




CP odd Higgs 


73.45 








75.25 






charged Higgs 


115.65 








116.92 






gluino 


725 








630 






charginos 


70 


260 






67 


232 




neutralinos 


52 


95 


129 


260 


48 


97 


113 232 


up squarks 


474 


614 


747 


771 


487 


603 


887 905 


down squarks 


510 


552 


749 


775 


513 


550 


893 909 


charged sleptons 


66 


338 


423 


448 


294 


550 


716 718 


sneutrinos 


327 




440 


440 


544 




713 714 
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Table 9: Model 4c - More Results at points IIIl and III2 (see points IIIl and III2 

on figures 4b and 4c ) 

Higgs, squark and slepton masses are treated in the same way as for the results quoted on 
Table VIII. 



Observable 


Predictions at point IIIl 


Predictions at point III2 


rriu/md \iGeV 


0.475 


0.432 




0.2337 








0.2344 








0.2325 








0.2331 






sin2Q; 


0.973 








0.977 






sin 2/5 


0.531 








0.516 






sin 7 


0.925 








0.918 






P 


-0.138 








-0.140 






'1 


0.331 








0.323 






tan/3 


55.39 








55.86 






CP even Higgses 


71.70 


112.36 






66.70 


112.05 




CP odd Higgs 


71.70 








66.70 






charged Higgs 


115.43 








112.22 






gluino 


447 








440 






charginos 


101 


218 






122 


274 




neutralinos 


62 


105 


174 


217 


67 


122 


251 272 


up squarks 


558 


677 


1443 


1460 


546 


666 


1440 1458 


down squarks 


597 


641 


1459 


1462 


574 


639 


1457 1460 


charged sleptons 


600 


1062 


1397 


1418 


564 


1051 


1397 1418 


sneutrinos 


1052 




1394 


1395 


1047 




1394 1395 
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Table 10: Model 4c - Varying the Lower Bound on 63 at point II Best fit values 
of selected GUT parameters and low energy observables are displayed for different lower 
bounds on €3. Note that the optimization procedure pushes 63 to its lower limit. The first 
line corresponds to the unrestricted case and the corresponding values are taken over from 
table 1^ for the sake of completeness. 







Mg X 10^^^ 


A 


as{Mz) 


X 10^ 


Mtop 


mb{Mb) 


Xtotal 


-4.89% 


24.362 


3.17 


0.807 


0.1151 


1.166 


175.7 


4.287 


0.731 


-3.00% 


24.628 


2.10 


0.808 


0.1154 


1.161 


176.1 


4.291 


0.959 


-2.00% 


24.767 


1.69 


0.802 


0.1156 


1.159 


176.1 


4.295 


1.279 


-1.00% 


24.907 


1.37 


0.813 


0.1158 


1.157 


176.6 


4.295 


1.694 


0.00% 


25.044 


1.11 


0.818 


0.1160 


1.155 


177.0 


4.296 


2.244 


+ 1.00% 


25.177 


0.89 


0.794 


0.1162 


1.153 


176.4 


4.304 


2.944 


+2.00% 


25.314 


0.73 


0.810 


0.1164 


1.151 


177.1 


4.302 


3.675 


+3.00% 


25.450 


0.59 


0.827 


0.1166 


1.149 


177.9 


4.303 


4.583 


+4.00% 


25.578 


0.48 


0.797 


0.1167 


1.147 


177.1 


4.309 


5.534 



Table 11: 

Model 4c - Sensitivity of the results of table ^ to the magnitude of the theoretical 
uncertainties In comparison to table |10|, the best fit values of the same GUT parameters 
and low energy observables are displayed for different lower bounds on €3, but this time for 
the uncertainty cr=0.1% for Mz, G^ and Oem, and the experimental value (T=130MeV for 
My/- (cr=0.5% for M^, My/ and Oem, and cr=1.0% for G^ was assumed in table [l^, as well 
as throughout this paper.) The first line corresponds to the case when €3 is unrestricted. 



£3 


1/"G 


Mg X 


A 


as{Mz) 


G^ X 10^ 


Mtop 




Xtotal 


-5.36% 


24.295 


3.53 


0.807 


0.1150 


1.166 


175.6 


4.286 


1.20 


-3.00% 


24.494 


2.55 


0.808 


0.1186 


1.165 


177.4 


4.418 


7.49 


-2.00% 


24.579 


2.21 


0.803 


0.1201 


1.165 


177.9 


4.471 


13.97 


0.00% 


24.759 


1.64 


0.818 


0.1229 


1.164 


179.6 


4.560 


34.02 
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FIGURE CAPTIONS 



Figure 1. 

Model 4 global analysis results in the tuq — M1/2 plane, for a fixed value of the Higgs pa- 
rameter 

a) fi{Mz) = 80 GeV , 

b) ^(Mz) =160 GeV , 

c) fj,{Mz) =240 GeV . 

Solid (dashed, double-dash-dotted) lines represent contour lines of constant = 15 (14,13) / 5dof. 
The star on figure la marks the point (mo=700GeV, Mi/2=240GeV). In the text it is referred 
to as point II and the corresponding results obtained at this point are listed in table ^ 

Figure 2. 

Results for Model 4 with one observable removed from the function (the corresponding 
standard deviation of that observable is inflated by a large factor), at fixed values of the 
Higgs parameter ij,{Mz)=80GeV and Mi/2=240GeV. Eight observables tried, (and their 
values obtained at the minimum of x^) ^-i'^ follows: as (0.1144) - dots, BR{b 57) 
(2.29- 10-2) - diamonds, Mb-M^ (2.84GeV) - solid squares, (1.50- lO'^) - open circles, 
Vub/Vcb (0.045) - triangles, mb(Mfe) (4.68GeV) - stars, Bk (1.12) - open squares, Kfe (0.0465) 
- solid circles. For comparison, the thick solid line represents the best values of Model 4 
with all observables present with their regular standard deviations as quoted in table |l|. No 
significant improvement of Model 4 is observed. Recall that leaving out an observable from 
the function means reducing the number of degrees of freedom by one. 

Figure 3. 

Model 4a global analysis results in the mo — M1/2 plane, for a fixed value of the Higgs 
parameter 

a) fi{Mz) = 80 GeV , 

b) fi{Mz) =160 GeV , 

c) i2{Mz) =240 GeV . 

Solid (dashed, double-dash-dotted) lines represent contour lines of constant = 6 (4,3) / 3dof. 
Figure 4. 

Model 4c global analysis results in the mo — M1/2 plane, for a fixed value of the Higgs 
parameter 

a) /i(M^) = 80 GeV , 

b) n{Mz) =160 GeV , 

c) ij{Mz) =240 GeV . 

Solid (double-dash-dotted, dotted) lines represent contour lines of constant = 6 (3,1) / 3dof. 
The stars mark the points (400,280), (700,240) and (1400,170). In the text, these are re- 
ferred to as points I, II, IIIl (point HI at /i=160GeV) and III2 (point HI at /i=240GeV). The 
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corresponding results and predictions obtained at these points are listed in tables § through 
i 

Figure 5. 

Stability of the predictions for sin2Q; and sm2(3, the parameters of the unitarity triangle, "x" 
symbols denote the points predicted by 

a) model 4c , 

b) model 4a 

at those regions in the SUSY parameter space where < 3 / 3dof . Dots in the figures 
represent the boundary of the allowed region resulting from a combined general Standard 
Model fit g for /s, = 200 ± AOMeV, Bb^ = 1.0, Bk = 0.75 ± 0.10. 



Figure 6a. 

The dependence of the quality of the fit on the GUT threshold correction €3, in model 4c. 
The thick solid line represents the results of model 4c with es unrestricted. The dotted line 
constrains €3 to be positive (and the optimization procedure always settles €3 very close to 
zero). Curves with solid squares, open circles and stars correspond to the results of the 
optimization performed with negative lower limit on -3%, -2% and -1%, respectively. 
Curves with open triangles, solid circles, open squares and solid triangles correspond to the 
cases with lower bound on €3 greater than zero, namely 1%, 2%, 3% and 4%, respectively. 
In all cases the optimization procedure tends to yield the values of 63 very close to the lower 
bound. Table [l^ displays the effect of different lower bounds on various GUT parameters 
and low energy observables at point II (i.e. at mo=700GeV ). 

Figure 6b. 

The same as in figure 6a, with a tighter theoretical uncertainty on precisely measured elec- 
troweak observables. In this figure the uncertainty for Mz, and is set to 0.1% and 
the actual experimental error 130MeV is assigned to a{Mw)] as opposed to the uncertainty 
of 0.5% for Mz, Mw and ttem, and cr(G'^)=1.0% , assumed in figure 6a. As before, in 
all restricted cases the optimization procedure tends to yield the values of e3 very close to 



the lower bound. Table 11 displays the effect of different €3 lower bounds on various GUT 



parameters and low energy observables at point II (i.e. at mo=700GeV ). 
Figure 7. 

The variation of the quality of the fit for different lower bounds on tree level pseudoscalar 
mass, in model 4c. The thick solid line represents the results of model 4c with m^*^^ >80GeV. 
When very close to this limit, the pseudoscalar mass receives negative corrections, and its 
one loop corrected value turns out to be at the experimental limit (see tables § and ^ for 
the results at point mo=700GeV). Curves with solid and open circles, solid squares, and 
stars correspond to the lower limit on m^*^^ of 130GeV, 200GeV, 300GeV and 500GeV, 
respectively. 
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Figure 4a 
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